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Abstract 

Diffractive production of D* ± (2010) mesons in deep inelastic scattering has been 
measured with the ZEUS detector at HERA using an integrated luminosity of 
44.3 pb _1 . Diffractive charm production is identified by the presence of a large 
rapidity gap in the final state of events in which a D* ± (2010) meson is recon- 
structed in the decay channel D* + — > (D° — > K~ii + )ii+ (+ charge conjugate). 
Differential cross sections when compared with theoretical predictions indicate 
the importance of gluons in such diffractive interactions. 
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1 Introduction 



Diffractive interactions in neutral current deep inelastic scattering (DIS) have been studied 
extensively at HERA [1-7]. Both inclusive and diffractive DIS cross sections rise more 
rapidly with energy than is the case for soft hadronic interactions [8], indicating the 
presence of a hard process to which perturbative Quantum Chromodynamics (pQCD) is 
applicable. Charm production is a key process for investigating the dynamics of diffractive 
DIS [9-14], since the charm-quark mass provides a hard scale and charm production is 
known to be sensitive to gluon-exchange processes in DIS [15,16]. 

Two contrasting approaches to describe diffractive DIS are considered in this paper: 

1. Resolved-Pomeron models: 

these are Regge-inspired, with an exchanged Pomeron having a partonic structure [17]. 
The evolution of the parton distributions of the Pomeron with Q 2 is described by 
the DGLAP equations [18]. Diffractive HERA data [6, 19] indicate a large gluonic 
component. In these models, charm is produced via boson-gluon fusion (BGF); 

2. Two-gluon-exchange models: 

these are based on decomposing the wave-function of the virtual photon in the proton 
rest frame into partonic Fock states — particularly qq and qqg [20]. These states then 
interact with the proton via colour-singlet exchange, the simplest form of which is 
the exchange of two gluons [21]. If the qq state dominates, charm production will be 
suppressed [22]. However, if the qqg state is dominant, similar production rates to 
those predicted by the resolved Pomeron model are expected [11, 13,23]. The qq and 
qqg configurations populate different regions of phase space. 

The "soft colour interaction" [24] and "semi-classical" [9] models give very similar pre- 
dictions for charm production to the two models described above and therefore are not 
considered separately. 

This paper describes the measurement of differential cross sections for D*^ production 
in diffractive DIS and the ratio with inclusive DIS D*^ production. The results are 
compared to resolved- Pomeron and two-gluon-exchange models. Similar results, including 
comparisons to the soft colour interaction, semi-classical and other models, have recently 
been published by the HI Collaboration [25]. 



2 Description of the experiment 

The integrated luminosity of 44.3 ± 0.7pb _1 used for this measurement was collected at 
the ep collider HERA with the ZEUS detector during 1995 - 1997, when HERA collided 
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27.5 GeV positrons with 820 GeV protons, giving a centre-of-mass energy of 300 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [26] . A brief outline of 
the components that are most relevant for this analysis is given below. Charged particles 
are tracked in the central tracking detector (CTD) [27], which operates in a magnetic field 
of 1.43 T provided by a thin superconducting coil. The CTD consists of 72 cylindrical drift 
chamber layers, organized in 9 superlayers covering the polar-angle[| region 15° < 9 < 164°. 
The transverse-momentum resolution for full-length tracks is <j(pt)/pt = 0.0058pT © 
0.0065 ©0.0014/pr, wither in GeV. The high-resolution uranium-scintillator calorimeter 
(CAL) [28] consists of three parts: the forward (FCAL), the barrel (BCAL) and the rear 
(RCAL) calorimeters. Each part is subdivided transversely into towers and longitudinally 
into one electromagnetic section (EMC) and either one (in RCAL) or two (in BCAL and 
FCAL) hadronic sections (HAC). The smallest subdivision of the calorimeter is called a 
cell. The CAL energy resolutions, as measured under test-beam conditions, are cr(E) / E = 
0.18/ VE for electrons and a(E)/E = 0.35/ \[E for hadrons (E in GeV). 

The position of positrons scattered with a small angle with respect to the positron beam 
direction was measured using the small-angle rear tracking detector (SRTD) [29], which 
also provides a means of correcting for any energy loss of the scattered positron due 
to the presence of inactive material. Complementing the SRTD is the rear presampler 
(RPRES) [30], which provides energy-loss information in RCAL regions outside the ac- 
ceptance of the SRTD. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > ejp, 
where the photon was measured with a lead-scintillator calorimeter [31] located at Z = 
-107 m. 



3 Kinematics of diffractive DIS 

The kinematics of the inclusive deep inelastic scattering of positrons and protons are 
specified by the positron-proton centre-of-mass energy, y^i, and any two of the following 
variables: Q 2 , the negative square of the four- momentum of the exchanged photon; y, 
the inelasticity; x, the Bjorken scaling variable and W, the centre-of-mass energy of the 
photon-proton system. 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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Additional variables are required to describe the diffractive process ep — > eXp, where 
X -> D^X'. These are: 

• t, the square of the four-momentum transfer at the proton vertex. Since t was not 
measured for the present data, all results discussed here are integrated over this vari- 
able; 

• x P = (M x + Q 2 )/(W 2 + Q 2 ), where Mx is the invariant mass of the hadronic final 
state, X, into which the virtual photon dissociates. This variable is the fraction of the 
proton's momentum carried by the exchanged colour- singlet system; 

• /3 = Q 2 / (M x + Q 2 ) can be interpreted within the resolved-Pomeron model as the 
fraction of the Pomeron's momentum carried by the struck parton. 

The above formulae for x P and /3 neglect the proton mass and assume t = 0. The variables 
are related by x = x F f3. 

To reconstruct the kinematic variables, both the final-state positron and the hadronic 
final state must be measured. The positron was identified using an algorithm based on 
a neural network [32]. The hadronic final state was reconstructed using combinations of 
calorimeter cells and CTD tracks to form energy-flow objects (EFOs) [4,33]. The DIS 
variables were reconstructed using the double angle (DA) method [34]. 

The mass of the diffractive system X was calculated from the EFOs as 

(?*)'- (?'-)'- (^)'- (?'-)'■ 

where the sum runs over all the EFOs in the event, excluding those associated with the 
scattered positron. 

4 Event Selection 

4.1 DIS selection and D* ± (2010) reconstruction 

The initial event sample was selected by identifying DIS events. The selection, both online 
and offline, was performed in the same manner as in the inclusive DIS D*^ study [15]. 
The only difference is in the kinematic region used, which for the present analysis is 
4 < Q 2 < 400 GeV 2 and 0.02 < y < 0.7. 

The D*^ selection cuts applied to reduce the combinatorial background differ somewhat 
from those used in the inclusive measurement, although the analysis is based on the same 
decay channel: D* + — > (D° — > K~7r + )irf (+ charge conjugate), where n s indicates the 
"slow" pion [35]. Reduction of the combinatorial background was achieved by requiring: 
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• the transverse momenta of any two oppositely charged tracks, assumed to be the K 
and the 7r from the D° decay, to be each greater than 0.5 GeV, and the transverse 
momentum of the slow pion from the D*^ decay to be greater than 0.12 GeV; 

• p(Ktt)/p(7v s ) > 8, where p(Kn) is the momentum of the candidate D° and p(ir s ) that 
of the slow pion. 

Since no particle identification was performed, the K and it masses were alternately at- 
tributed to the decay products of the candidate D° meson. Only D° candidates that had 
an invariant mass between 1.80 GeV and 1.92 GeV were subject to the mass difference re- 
quirement 0.143 < AM < 0.148 GeV (AM = M(Kmr s ) - M(Ktt)). Further requirements 
placed on the D** candidate were 1.5 < pt^D*^) < 8.0 GeV and ^(Z?**)] < 1.5. 

After applying these requirements, a signal of 1720 ± 63 D*^ mesons was obtained, using 
the fit procedure described in Section |4~3j 

4.2 Diffractive selection 

A key characteristic of a diffractive event is the presence of a large rapidity gap between 
the scattered proton, which remains in the forward beampipe, and the hadronic system 
X from the dissociated virtual photon. Such events were selected using a cut on ?7 max , 
the pseudorapidity of the most-forward EFO with energy greater than 400 MeV in the 
event [36]. Figure [l|(a) shows the distribution of r) max for events containing a D* ± candi- 
date. The non-diffractive events deposit energy around the FCAL beamhole, the edge of 
which is located at a pseudorapidity of about four units, producing the peak at ?7 max ~ 3.5. 
Below that value, the contribution from non-diffractive interactions falls exponentially [4], 
leaving a plateau at lower 7/ max values, which is the signature of diffractive interactions. 
The shaded area in Fig. 0(a), obtained from the non-diffractive Monte Carlo simulation 
described in Section |5], exhibits an exponential fall-off at low r] max and has no events with 
?7max < 2, confirming that the plateau observed in the data corresponds to diffractive 
events. A cut of ?7 max < 2, i.e. a gap of at least two units of pseudorapidity, was used to 
select diffractive events. 

There are two main implications of the use of the ^ max selection. First, this method defines 
the rapidity gap with respect to the forward edge of the calorimeter and therefore cannot 
distinguish diffractive events in which the proton remains intact from those in which the 
proton dissociates into a low-mass system, whose decay products remain in the forward 
beampipe. To correct for this, a proton-dissociative contribution of (31 ± 15)% [4] has 
been subtracted from all measured cross sections. Secondly, the r/ max method restricts the 
range of x P values; to account for this, a cut of x P < 0.016 was applied. 
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The presence of charm in the diffractive events automatically sets a lower limit on Mx- 
This, in turn, places an upper limit on the value of j3 of f3 max ~ 0.96. Since, in addition, 
the acceptance falls steeply at high (5, only events that satisfy (3 < 0.8 were retained for 
further analysis. 

4.3 D*± Signal 

The AM distribution after all selection cuts is shown in Fig. |l](b). A clear D*^ signal is 
evident over a small background. 

To determine the number of D*^ candidates, the AM distribution was fitted using a 
Gaussian function for the signal and the form 

a (AM - m w ) b e c{AM - m ^ 

for the background, where a, b and c, as well as the mean, width and normalisation of 
the Gaussian function were free parameters of the fit. The number of D*^ candidates 
resulting from an unbinned log-likelihood fit was 84 ± 13. The mean value of AM for the 
peak was 145.67 ± 0.14 MeV, consistent with the Particle Data Group (PDG) value [37]. 
The width of the Gaussian was 0.93 ±0.16 MeV, in agreement with the detector resolution. 
The fit is good, as seen in Fig. 0(b). The same fit procedure was then followed in each 
cross section bin to obtain the number of D*^ 1 candidates and gave satisfactory results in 
each case [38]. The systematic uncertainties relating to the extraction of the number of 
candidates are discussed in Sect. |[ 

The Kir invariant-mass distribution was fitted with a Gaussian function for the signal 
and a simple polynomial to describe the background. The fit to the M(Kir) distribution 
yielded a D° mass of 1865.7 ±2.1 MeV and a resolution of 15 ± 2 MeV. The former is 
in good agreement with the PDG value [37], while the latter is in agreement with the 
detector resolution. The fit results for both the AM and the M(Ktt) invariant-mass 
distributions are comparable to those found in the inclusive dataset [15,38]. 

5 Monte Carlo simulation and reweighting 

A GEANT-based [39] Monte Carlo (MC) simulation was used to calculate selection effi- 
ciencies and correction factors. Two different diffractive event generators were used: RIDI 
v2.0 [40] for evaluating the nominal correction factors and RAPGAP v2. 08/01 [41] as a 
systematic check. For all the MC samples, events with at least one D*^ decaying in the 
appropriate decay channel were selected and passed through the standard ZEUS detector 
and trigger simulations as well as the event reconstruction package. 
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The RIDI generator is based on the two-gluon-exchange model developed by Ryskin [40]. 
All stages of the parton fragmentation and hadronisation were simulated with the Lund 
string model [42], as implemented in JETSET [43]. The simulation included QED radia- 
tive corrections. Separate samples of 7* — > qq and 7* — > qqg interactions were generated. 
The CTEQ4LQ [44] proton parton distribution functions were used as input and the 
charm mass (m c ) was set to 1.35 GeV. The acceptance corrections have little sensitivity 
to the value chosen for m c . 



The distributions of the variables Q 2 , x P and j3 from the qq and qqg RIDI MC samples are 
compared to the data in Fig. ^. Neither MC sample alone reproduces the measurements, 
particularly the f3 distribution. Therefore, the two RIDI MC samples were combined by 
fitting the fraction of qq events in the three bins of Q 2 used to extract the differential 
cross sections. The x P distribution in each Q 2 bin was used in making the fit. It was 
found that neither reweighting in (3 nor in x P separately gave a good description of the 
detector-level distributions in the data. The fraction of qq events increases with Q 2 from 
approximately 4% in the lowest-Q 2 bin to around 50% in the highest-Q 2 bin. This mixed 
sample, shown as the hatched histograms in Fig. ^ is in satisfactory agreement with the 
data. The acceptance in the kinematic region used for the cross-section measurement was 
17.3%. 

The RAPGAP simulation provides a rather general framework for the generation of 
diffractive events. The resolved-Pomeron option was used, in which charm quarks are 
produced via the leading-order BGF process. Charm fragmentation was carried out using 
the Peterson fragmentation model with the parameter e c set to 0.035 [45]. The fragmenta- 
tion process was simulated using the colour-dipole model implemented in ARIADNE [46] 
and the hadronisation was carried out according to the Lund string model. The sample 
was generated assuming a gluon-dominated Pomeron, with a gluon distribution peaked 
close to /3 — 1 [5,6]. The Pomeron intercept was set to a P (0) = 1.20 and the charm 
mass was set to 1.35 GeV. The QED radiative corrections are not available for diffractive 
charm production in RAPGAP. 

Similar discrepancies to those observed in the original RIDI samples were seen when 
comparing the RAPGAP MC sample to the data, with the exception that the x P distri- 
bution was well described. To correct for the discrepancies, the RAPGAP sample was 
simultaneously reweighted in logQ 2 and Mx, after which agreement was obtained in all 
distributions. 

The RAPGAP generator, which also simulates non-diffractive interactions, was used to 
produce a non-diffractive D*^ sample for the extraction of the ratio of diffractive to 
inclusive Z)* ± production (see Section |7~2]) . The parameters used were the same as those 



used in the measurement of the inclusive DIS D*^ 1 cross sections [15]. 
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6 Systematic uncertainties on the cross section mea- 
surements 



The major sources of systematic uncertainties and their effect on the measurement of the 
cross section are: 

• the selection of inclusive DIS events. These systematic uncertainties were derived in 
the same way as for the inclusive DIS D*^ analysis [15] and resulted in an overall 
variation in the cross section of lg%; 

• the selection of D*^ candidates. The minimum transverse momentum of tracks used 
in the D* ± reconstruction was raised and lowered by 15% and the threshold on the 
momentum-ratio p(Ktt) /p(tt s ) was increased by 0.5 unitsQ These variations yielded a 
combined uncertainty of ±6%; 

• the selection of diffractive events. The 77 max requirement was varied by ±0.2 units 
and the EFO energy threshold applied in the reconstruction of 77 max was varied by 
±100 MeV. The combined effect of these changes was lg%; 

• the model dependence of the corrections. This uncertainty was estimated using the 
reweighted diffractive RAP GAP sample instead of the mixed RIDI sample. This 
change resulted in a variation of —9%; 

• the reweighting procedure. The qq fraction in each Q 2 bin was separately varied up 
and down by its uncertainty as determined from the fit (see Section ||). The resulting 
deviations were summed in quadrature, giving a variation of tfflo- 

These systematic uncertainties were added in quadrature separately for the positive and 
negative deviations from the nominal values of the cross section to determine the overall 
systematic uncertainty of i 16 %. These estimates were also done in each bin in which 
the cross section was measured. The sources of systematic uncertainty relating to the 
extraction of the number of D* ± candidates were also considered, but were found to be 
negligible in comparison to those listed above. 

The overall normalisation uncertainties arising from the uncertainties on the luminosity 
measurement and the D*^ and D° branching fractions were not included in the sys- 
tematic uncertainty. The uncertainty arising from the subtraction of proton-dissociative 
background, quoted separately, is ±22% [4]. 

7 Results 



2 A trigger requirement did not permit the decrease of this momentum-ratio cut. 
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7.1 Cross sections 



The differential cross sections in any given variable £ were determined using 

da _ N D (1- /p diss ) 

d£ ~ A ■ £ ■ B ■ 5i ' U 

where Np is the number of D*^ candidates fit in a bin of width <5£, A is the correction 
factor (accounting for acceptance, migrations, efficiencies and radiative effects) for that 
bin, £ is the integrated luminosity and B = 2.59 ± 0.06% is the total D* + — > D tt^ — > 
(if~7r + )7r+ branching ratio [37]. The quantity /pdiss is the fraction of proton-dissociative 
background. All other sources of background were neglected. 

Figure |3| shows the measured differential D*^ cross sections with respect to the kinematic 
variables Q 2 , W, x P , j3, px{D* ± ) and ?7(_D* ± ). The values are also given in Table |l[ The 
data exhibit a sharp fall-off as a function of Q 2 and priD^). The events are concentrated 
at low j3, but are uniformly distributed in x P , within the large uncertainties. The W 
dependence is mainly determined by the pr(-D* ± ) and ?7(D* ± ) restrictions. 

The cross section for diffractive D*^ production in the kinematic region 4 < Q 2 < 

400 GeV 2 , 0.02 < y < 0.7, 1.5 < M- **) < 8 GeV - W- **)! < L5 > x *> < °- 016 and 
(3 < 0.8 is 

(?ep^eD*±x> P = 291 ± 44(stat.)l^(syst.) ± 63(prot. diss.) pb, (2) 

where the last uncertainty arises from the subtraction of the background from proton 
dissociation. 



7.2 Ratio of diffractive to inclusive D production. 

The ratio of diffractive to inclusive D*^ production was measured for x < 0.0128. This 
limit is equivalent to the x P and f3 requirements imposed on the diffractive sample. In the 
kinematic region 4 < Q 2 < 400 GeV 2 , 0.02 < y < 0.7, 1.5 < priP**) < 8 GeV, |?7(-D* ± )| < 
1.5 and x < 0.0128, the inclusive DIS D*^ 1 cross section is a ep ^ eD *± Y = 4.83 ± 0.18(stat.) nb, 
where Y is the complete hadronic final state except for the tagged D*^ 1 meson. This value 
is consistent with an earlier ZEUS measurement [15], taking into account the differences 
between the two kinematic regions. The ratio of diffractive to inclusive DIS D*^ 1 produc- 
tion is then defined by 

= a e p-,eD*±x'p(xjp < 0-016, g < 0.8) 
° ~ a ep ^ eD *± Y (x < 0.0128) ' 1 ' 

It is assumed that all the systematic uncertainties, except those relating to the diffractive 
selection and the MC reweighting, cancel. This is a reasonable approximation given that 
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the cross section results in the diffractive case are statistically limited. The measured 
ratio is therefore 

R D = 6.0 ± 0.9(stat.)+^(syst.) ± 1.3(prot. diss.)%. (4) 
Figure HI shows Rd clS 8b function of Q 2 and W. 

8 Discussion 

Three models are compared to the measured cross sections: (1) the resolved-Pomeron 
model, as implemented in the fits to HERA data made by Alvero et al. (ACTW) [47], 
(2) the two-gluon-exchange "saturation" model of Golec-Biernat and Wiisthoff [48], as 
implemented in the SATRAP MC generator and interfaced to RAPGAP [49] and (3) the 
two-gluon-exchange model of Bartels et al. (BJLW) [14,50], as implemented in RAPGAP. 

The version of SATRAP used is the same as that used in the recent ZEUS publica- 
tion on the study of hadronic final states in diffraction [7], where it is referred to as 
"SATRAP-CDM" . Although the SATRAP and BJLW predictions are both based on two- 
gluon exchange, they differ in the treatment of the qqg final state, which is an important 
contributor to charm production. In the SATRAP-CDM model, only configurations in 
which the transverse momenta satisfy kr(q), k^q) > kr(g) are included in the calculation 
of the qqg final state, while all configurations with kx{g) > 1 GeV are used in the BJLW 
model. 

The predicted cross sections from three of the fits by Alvero et al. are shown in Table |^. 
The calculations were made assuming m c = 1.45 GeV and the Peterson fragmentation 
model (with the parameter e c = 0.035) for the charm decay. The probability for charm to 
fragment into a D*^ meson was taken as 0.235 [51] and the renormalisation and factori- 
sation scales were set to /i# = fip = \/Q 2 + 4m 2 .. Varying the charm mass by ±0.15 GeV 
gave a t\*% 

variation in each cross section. The theoretical uncertainties relating to the 
scale and the value of the Pomeron intercept are comparable to that arising from the 
charm mass. Fits B and D, which assume a gluon-dominated Pomeron, are favoured by 
the data. The fit SG, which assumes a Pomeron dominated by gluons with a "superhard" 
density, predicts too small a cross section. A quark-dominated Pomeron (not shown) 
predicts a cross section too small by two orders of magnitude. 

Table |2| also shows the predicted cross sections from SATRAP-CDM and the BJLW model, 
along with the qq contribution from the BJLW calculations. The same assumptions are 
made in these calculations as were made in the ACTW case, except that here (1r — Hf = 
Q. The BJLW qq contribution is clearly too small, while the sum of the qq and qqg BJLW 
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contributions gives good agreement with the data. The SATRAP-CDM and ACTW fit B 
predictions are smaller than the data, but not significantly so. 

Figure | compares the predictions of ACTW fit B, SATRAP-CDM and BJLW, along 
with the qq BJLW contribution, with the measured differential D*^ cross sections. The 
ACTW fit B and fit D predict similar shapes for the differential distributions and therefore 
no comparisons with fit D are made. The ACTW and SATRAP-CDM predictions are 
in reasonable agreement with the measured differential distributions, with the exception 
of the (3 distribution, where both predictions undershoot the data at high (3. The qq 
contribution from the BJLW calculations clearly fails to describe the shape of the distri- 
butions, particularly that for (3. The sum of the two BJLW contributions gives a good 
representation of all the measured distributions. 

The ratio, Rd, agrees with that determined in inclusive diffraction [4], indicating that 
charm production is not suppressed in diffractive DIS, contrary to the expectations of 
some early models [22]. Figure |] shows that Rd is consistent with being independent of 
Q 2 and W. 

The HI Collaboration recently published similar results [25]. The kinematic region in 
that study is significantly different to that used here, making it difficult to compare the 
results directly. However, comparisons with Monte Carlo models indicate that the two 
sets of results are consistent. 

9 Summary 

The cross section for diffractive D*^ production in the kinematic region, 4 < Q 2 < 
400 GeV 2 , 0.02 < y < 0.7, x P < 0.016, (3 < 0.8, 1.5 < M- **) < 8 GeV and 1^(^)1 < 
1.5 has been measured to be 291 ± 44(stat.)^47(syst.) ± 63(prot. diss.)pb. Differential 
cross sections have been compared to the predictions of different models of diffractive 
charm production. The resolved-Pomeron model of Alvero et al. is below the data at 
high (3. In the case of the two-gluon-exchange models, the qq contribution alone gives too 
small a cross section. Inclusion of the qqg contribution results in a good description of 
the data by the BJLW prediction, while the SATRAP-CDM prediction undershoots the 
data, particularly at high (3. 

The ratio of diffractive D*^ production to inclusive DIS D* ± production is Rd = 6.0 ± 
0.9(stat.)lS:?(syst.)±1.3(prot. diss.) %. This result is in agreement with the correspond- 
ing ratio for inclusive diffraction, indicating that charm production in diffraction is not 
suppressed with respect to light-flavour production. The ratio Rd is independent of Q 2 
and W within the uncertainties. 
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da/dQ 2 


Q 2 Range ( GeV 2 ) 


Q 2 ( GeV 2 ) 


a (pb) 


4 - 12 


6.7 


131 ± 32±|| ± 28 


12 - 25 


17.8 


95 ±21^ ±21 


25 - 400 


67.7 


64 ± 18^° ± 14 


da/dW 


W Range ( GeV) 


W ( GeV) 


tr(pb) 


50 - 130 


97.4 


121 ±41+^ ±26 


130 - 170 


149.3 


84 ± 20+^ ± 18 


170 - 250 


213.9 


97 ±27^ ±21 


da/dx P 


x P Range 


Xp 


*(pb) 


- 0.0045 


0.0031 


99 ± 26i^ ± 22 


0.0045 - 0.009 


0.0066 


68 ± 17+^ ± 15 


0.009 - 0.016 


0.0123 


150 ± 34±?gg ± 33 


da/d(3 


P Range 


P 


a (pb) 


- 0.1 


0.043 


129 ± 31+^ ± 28 


0.1 - 0.3 


0.19 


72 ± 16^ 6 ± 16 


0.3 - 0.8 


0.50 


81 ± 20l? 2 ± 18 


da/dpriD^) 


M- **) Range (GeV) 


MD**) (GeV) 


a (pb) 


1.5 - 2.4 


1.9 


149 ± 491 J| 2 ± 32 


2.4 - 3.6 


2.8 


81 ± 20+?3 ± 18 


3.6 - 8.0 


4.4 


52 ± lllf ± 11 


d(x/dr/(D* ± ) 


r7(D* ± ) Range 




a (pb) 


-1.5 - -0.65 


-0.98 


106 ± 28±|o ± 23 


-0.65 - 0.1 


-0.24 


95 ±21+^ ±21 


0.1 - 1.5 


0.79 


76 ± 23j$ ± 17 



Table 1: Values of the differential cross sections with respect to Q 2 , W, x P , P, 
Pt{D*^) andri{D* ± ). The following quantities are given: the range of the measure- 
ment; the value at which the cross section is quoted and the measured cross section. 
The first uncertainty quoted represents the statistical uncertainty, the second the 
systematic and the third the uncertainty arising from the subtraction of the proton 
dissociation background. 
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Cross Section (pb) 


Data 


291 ± 44+^ ± g3 


ACTW FIT B 


187 


ACTW FIT D 


401 


ACTW FIT SG 


87 


SATRAP 


185 


BJLW qq 


79 


BJLW qq + qqg 


297 



Table 2: Comparison of the measured cross section with those predicted by the 
models described in the text. No errors are quoted given the considerable flexibility 
in the predicted cross sections depending upon the choice of input parameters for 
the models. 
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Figure 1: (a) The ?7 max distribution for DIS events with a D*^ candidate. The 
shaded histogram is the prediction from the RAPGAP non-diffractive Monte Carlo 
simulation (see section 5). (b) The AM distribution for events with r/ max < 2. Only 
the combinations whose values of M (Kir) lie in the signal region are included. The 
curve shows the unbinned log-likelihood fit described in the text. 
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Figure 2: The Q 2 , x F and [3 distributions at the detector level for the data and 
the RIDI MC samples. The data sample is shown as the solid points and the qq and 
qqg RIDI samples are shown separately as the solid and dashed histograms. The 
mixed sample, produced according to the procedure described in the text, is shown 
as the hatched histograms. All the MC distributions have been normalised to the 
number of events in the data. 
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Figure 3: Differential cross sections for D*^ production in the kinematic region 
described in the text. The cross sections are shown as a function of (a) Q 2 , (b) 
W, (c) x P , (d) (3, (e) pt^D^) and (f) ^(D*^). The inner bars show the statistical 
uncertainties, while the outer bars indicate the statistical and systematic uncertain- 
ties added in quadrature. The shaded bands show the uncertainty coming from the 
subtraction of the proton-dissociation background. The histograms correspond to 
the different models described in the text. 
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Figure 4: The measured ratio of diffractively produced D*^ mesons to inclusive 
D*^ meson production, R^, as a function of Q 2 and W. The inner bars indicate 
the statistical uncertainties, while the outer ones indicate the statistical and system- 
atic uncertainties added in quadrature. The shaded band indicates the uncertainty 
arising from the subtraction of the proton- dissociation background. 
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